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ABSTRACT 

Jets are ubiquitous in astronomy. It has been conjectured that the existence of 
jets is intimately connected with the spin of the central object and the viscous an- 
gular momentum transport of the inner disk. Bipolar jet-like structures propelled 
by the viscous torque on a spinning central object are also known in a completely 
different context, namely the flow in the laboratory of a viscoelastic fluid. On the 
basis of an analogy of the tangled magnetic field lines of magnet ohydro dynamic 
(MHD) turbulence to the tangled polymers of viscoelastic polymer solutions, we 
propose a viscoelastic description of the dynamics of highly turbulent conductive 
fluid. We argue that the same mechanism that forms jets in viscoelastic fluids 
in the laboratory may be responsible for collimating and powering astrophysical 
jets by the angular momentum of the central object. 

Subject headings: MHD — turbulence — stars: winds, outflows — galaxies: jets 

1. Introduction 

Viscoelastic behavior is characteristic of fluids such as polymer suspensions in which long 
molecular chains become entangled. The deformation and disruption of this entanglement 
that occurs when the fluid is subjected to shear, and the relaxation to a new entangled state, 
give such fluids a "memory" that ordinary Newtonian fluids do not have. Egg whites are 
one common household example. These fluids are in a sense midway between being elastic 
solids and Newtonian fluids, being more one or the other depending upon the ratio of the 
relaxation time of the entangled polymers to the deformation time scale. This ratio is known 
as the Weissenberg number We. Viscoelastic behavior can cause flows that are qualitatively 
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quite different from the Newtonian case, such as the Weissenberg effect, in which a fluid 
climbs a spinning rod. 

We propose that the effects of "mesoscale" magnetic fields in some instances of MHD 
turbulence can be modeled by additional viscous terms, and that the structure of these terms 
is similar to that of a viscoelastic fluid. (It has come to our attention that this has also been 
conjectured in a recent paper by Ogilvie (Ogilvie 2001), who demonstrated that a viscoelastic 
model of MHD turbulence can dramatically affect the dynamics of eccentric accretion disks.) 
It has long been noted that the dynamics of vortex filaments of fully-developed turbulence 
bear a superficial likeness to a network of polymers; see Chorin (1988). There is also an 
obvious analogy between the tangled polymers of a polymer fluid and the tangled field lines 
in MHD turbulence in the case that most of the field energy is on intermediate to small scales. 
However, in contrast with polymers which are typically conserved, the magnetic energy on 
small scales is not conserved and may dissipate. We therefore restrict our attention to flows in 
which the rate of shear changes slowly with respect to the dissipation time scale. We assume 
that the large-scale {i.e. mean-field) magnetic field is rather weak, but the magnetic field on 
intermediate scales is advected by hydrodynamic turbulence (Ott 1998). The particular case 
we consider is turbulence driven by shear, such as in a Balbus-Hawley magnetic shearing 
instability (Balbus & Hawley 1998) in the case that the net magnetic flux is negligible. 

One of the more startling effects seen in viscoelastic fluids is the reversal of the secondary 
flow around a rotating sphere at moderate rotation speeds. A Newtonian fluid will flow 
inwards along the rotation axis and be expelled in the equatorial plane, but in a viscoelastic 
fluid, the secondary flow is exactly the opposite. We propose this phenomenon as a model 
for the production of astrophysical jets powered in part by the angular momentum of the 
central object. 



In MHD, the magnetic field evolves according to the vector advection-diffusion equation 



where B is the magnetic field, v is the fluid velocity and r] is the magnetic diffusivity, which 
is determined by the conductivity of the fluid. If the diffusivity rj is zero, the flux-freezing is 
perfect. The contribution of magnetic fields to the stress tensor is given by 



2. MHD "rheology 



dtB = B ■ Vv — v • VB — B(V • v) + ??V 2 B, 
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Let us average spatially over some "subgrid" length scale, and concentrate on the curvature 
term. Under perfect flux- freezing, the tensor = {BiBj) can be shown from eq. (1) to 
evolve according to 



where Lij = diVj and D t is the usual advective derivative. The first three terms in eq. (3) 
are known collectively in rheology as the upper-convected invariant derivative T>^ of the 
tensor M^. Insofar as the viscous stress in a non-Newtonian fluid does not depend only on 
the velocity gradient L^, but on its history as well as higher-order derivatives of the velocity, 
the concept of the advective derivative d t + v • V must be extended to respect the material 
frame indifference of the stress tensor. The upper-convected invariant derivative £>j u ' of the 
viscous stress tensor is a common material derivative in rheology, and it is used in many 
rheological models. 

According to flux- freezing, as expressed in eq. (3), on- and off-diagonal components of 
the magnetic stress in a highly conducting plasma threaded with a small-scale magnetic field 
will grow in response to large-scale shear. A simple analogy is a square piece of rubber on 
which we have drawn a tangled, statistically homogeneous and isotropic scribble. If we pull 
the upper edge of this sheet to the right, say, and the lower edge of the sheet to the left, the 
area of the square will remain the same, but the scribble is now no longer isotropic. There 
is an overall tendency of the scribble to be composed of lines oriented diagonally from the 
lower left corner to the upper right corner. 

With this picture in mind, suppose that at some time t — the individual components 
of the magnetic field are distributed as Gaussian random variables with equal variance a 
and zero mean. In a Cartesian reference frame, given a uniform shear 7 in the x direction, 
with perfect flux-freezing, the mean field will evolve according to 



The field lines are in reality constantly shifting and being arranged on small scales by small 
turbulent eddies, which we characterize by a relaxation time s. This time incorporates 
the effects of the turbulent stretching and twisting that generates the field, as well as the 
dissipative processes that are governed by the magnetic diffusivity 77. In the case we consider, 
the magnetic field draws its energy from the shear, so it is natural to assume that 1/s is 
comparable to the rate of shear 7. We therefore expect that the product 7s, which is the 
Weissenberg number (We), is of the order of magnitude of one or greater. (Ogilvie finds 
We ~ 1 based on a comparison with simulations.) If the turbulence is driven by other 
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processes, such as convection, the effective Weissenberg number might be very different from 
1. 

In the case of steady shear, we approximate the form of the tensor (BiBj) by replacing 
t in the above expression with the relaxation time s. The off-diagonal stress is proportional 
to 7s = (We), with the constant of proportionality determined by the viscosity prescription. 
There is a normal stress difference in the direction of shear t xx — r yy = (We)r xy . This 
form of the stress tensor, modulo the magnetic pressure, is very similar to the second-order 
approximation to a viscoelastic Weissenberg fluid, with zero second normal stress difference 
Tyy — t zz . The non-Newtonian viscous stress in such a fluid may be written in terms of 
the Rivlin-Erikson tensors and as r = /j,A^ + a\A^ + a 2 (A^) 2 where is the 
usual viscosity and the second-order coefficients «i and a 2 satisfy a 2 = —2a\. The first 
Rivlin-Erikson tensor A^ is simply the symmetrized velocity gradient, Afj = Lij + Lji, 
and subsequent tensors A^ obey the recursion relation A!™" 1 " 1 ' = D t A^ + A^L kj + A^L ki . 
The nth Rivlin-Erikson tensor gives the nth time derivative of the length of an infinitesimal 
material line |dx| embedded in the fluid, so that d^ |<ix| 2 = A^dxidxj. A second-order 
fluid is thus simply an approximation to a non-Newtonian fluid that is second-order in the 
velocity gradient, having contributions to the viscous stress that are proportional to the 
second spatial derivative of the velocity as well as the square of the first spatial derivative of 
the velocity, in addition to the normal Newtonian viscous stress which is proportional to the 
(symmetrized) first spatial derivative of the velocity. As has been noted by Bohme (1987), 
the term "second-order fluid" is thus somewhat unfortunate since it is the flow, and not the 
fluid, which is really second-order. If the material derivative of the rate-of-strain tensor itself 
changes on a time scale comparable or faster than s then a more sophisticated model must 
be brought to bear. 

In a cylindrical coordinate system with azimuthal shear, such as in the fluid sur- 
rounding a rotating sphere, the normal stress difference contributes to a global hoop-stress 
T ee , which acts like a rubber-band, creating a radially-inward "non-Newtonian" specific 
force / nn = —Tee/{rp). In the laboratory, the hoop-stress-induced forces cause a radially- 
inward flow close to the equatorial plane (Giesekus 1963; Geisekus 1994). This has also 
been shown mathematically up to various orders in We and Re for a second-order fluid 
(Giesekus 1963; Bohme 1987; Joseph 1990; Geisekus 1994) as well as to leading order for 
a viscoelastic Oldroyd model B' fluid (Thomas & Walters 1964) which has a spectrum of 
relaxation times. This inflow pushes and collimates the intervening material into bipolar 
axial outflows for moderate We and Reynolds number Re. 
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3. 



Jet model 



The astrophysical analog we propose is a spinning massive central object immersed in 
a viscous thick disk with a tangled magnetic field. The solutions in the rheology literature 
do not map directly onto astrophysical solutions because gravitational, compressibility, and 
radiative effects are important in the latter. Furthermore, the viscosity in the astrophysical 
case is probably not constant as we assume here. We therefore confine ourselves to general 
dimensional arguments. We will allow the central object to have an extended magnetosphere 
that acts like a rigid rotator, so that the effective radius of the object Rq may in fact 
be much larger than the actual radius of the object, depending upon the strength of the 
magnetosphere of the object. We assume that the disk and central object are viscously 
coupled so that Qq = Q(Rq), and we define Q to be the angular speed of the central object 
Q normalized to the Keplerian frequency at the object's effective surface Qk(Ro)- Note 
that we may have u> > 1. The central object rotates in the same direction as the disk. Also 
we define r = R/R . Define r m to be the dimensionless radius of the thick disk, and let the 
angular frequency of the disk f2 scale as R~ q . Our conclusions are not sensitive to this last 
assumption, but it provides a nice framework to discuss the viscoelastic effects. Note that if 
we take g>3/2foru;>l and q < 3/2 for Q < 1, the thick disk will match to a Keplerian 
disk at r 9 ~ 3/2 = uj 

The structure of thick disks is sensitive to the viscosity prescription, and the "correct" 
viscosity is unknown. We assume that Rq < H < r m Ro where H is the thickness of the thick 
disk. As an upper and lower bound, we set 



where (3 is between 1 and and a ~ 1 as usual. The traditional parameter a is thus 
redundant here. The effective Reynolds number as determined by Rq and Qq is then Re = 
l/(qa(3), which is on the order of 1 or less. 

The viscous stress is r r g = pu(RQ'). The specific inward "non-Newtonian" radial force 
f nn due to the normal stress difference is then 



It must not be forgotten that the viscous and hoop stresses occur in proportion to the 
magnetic pressure. From eqs. (2) and (4), and recalling that t t q = v(RQ!), we find 
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The specific centrifugal inertial force is /i ner t = & 2 R and the specific gravitational force 
is /grav — —^ij^R, so that the inward hoop stresses will dominate inertial forces for r 2 < 
(We)q 2 a[3, and they will dominate gravitational forces for r^ 29-1 ) < (We) q 2 a (3 to 2 . However, 
it is important to note that, in the incompressible case, or for a polytrope, the specific forces 
/grav , /inert; and / pre ss are all derivable from potentials, whereas / nn is not. Note that the 
specific pressure force /press = p _1 VP may be written as the gradient of some function so 
long as there is a one-to-one relation between P and p. We expect that there is typically 
no realistic combination of gravitational, centrifugal, or pressure forces that will balance the 
non- Newtonian hoop stresses for an arbitrary choice of vq(R, z), regardless of the value of (3. 
For the choice of d z v e = 0, such as in a thin disk, the hoop stress / nn depends only on R and 
thus may be written / nn = V<& for some $, but for a thick disk this will not typically be the 
case. 

In contrast with azimuthal viscous forces which dissipate kinetic energy, f nn -v is positive 
for a radial inflow, so that the hoop stresses do positive work on the fluid. There are no 
hoop stresses along the rotation axis of the system, so the energy gained by the fluid is free 
to escape in this direction. Lastly, as in the laboratory case, it is hypothesized that these 
hoop stresses naturally confine and collimate the flow, if the thick disk is sufficiently larger 
than the object radius Rq. 

The hoop stresses are powered in part by the mechanical accretion luminosity, and the 
specific work performed by them for a fluid element in the equatorial plane that moves from 
the outer radius R m to i? is 

rRm i CM 

W= f nn dR= qa /3(We) — (l - r m 2 «) (8) 

JR 1 -tto 

For (5 >> 1 this work can be much larger than the binding energy; some of this energy will 
be viscously dissipated, and some is work against the pressure gradient. 

Extrapolating from the laboratory case for moderate We and Re, the launch speed of 
the jet is expected to be of the order of 

where C2 is a constant of order unity. In escaping to infinity, the jet will lose its binding 
energy GM/R Q = (RqQo/oj) 2 , but it will be propelled by the magnetic pressure P mag / p given 
in eq. (7). If the thick disk is substantially thicker than the object effecive radius Rq then the 
jet may also suffer some viscous drag, which we estimate by Poiseuille's law. In traversing 
the thickness of the disk if a jet of radius aR (where we expect a ~ 1) will lose an amount 



(9) 



7 



of energy per unit mass of the order of vv- ]et H/ (aRo) 2 = C^HRqVIq where C 3 is of order unity. 
We thus very roughly estimate the jet speed at infinity as 



The speed of the jet naturally scales as Rq£Iq, and increases with the speed of rotation u>, 
which must be close to breakup for /3 ~ 1 for the jet to reach infinity. We caution however 
that further work is required to find to what degree the above admittedly very approximate 
relationship holds in the astrophysical case. This expression is intended only to show roughly 
the relative importance of the various effects we have considered. We have not considered 
any relativistic effects, and we have not addressed the importance of the accretion efficiency, 
nor any effects associated with the central object having an event horizon. Accretion onto a 
black hole will certainly complicate matters very much. As it stands, therefore, we believe 
our present arguments apply perhaps best to the case of jets from protostellar objects. 



We have proposed that the fluid dynamics of a highly-conductive turbulent medium 
with tangled magentic fields is similar in some respects to the viscoelastic fluid dynamics 
of a suspension of tangled polymers. We discuss the production of jet-like structures in 
viscoelastic fluids induced by a spinning sphere, and outline an analogous model for the pro- 
duction of astrophysical jets propelled by a spinning central object. The processes outlined 
above are thus hypothesized to propel jets by tangled magnetic fields, as has been previously 
suggested by Heinz & Begelman (2000). 

We believe that this tentative model for the acceleration of jets has some advantages 
over previous models. Firstly, the jet process does not rely upon the insertion "by hand" of 
any large-scale magnetic field. The jet also appears to be self-collimating, and it is a solid- 
core jet, rather than a hollow-core jet. As pointed out by Li (2001), tangled magnetic fields 
may also stabilize the jet from the long-wave mode kink instability that plagues standard 
MHD jet models. Lastly, the jet mechanism we describe is powered in part by the angular 
momentum of the central object. We hope to improve upon this work with more detailed 
calculations in a future paper. 

The majority of this work was completed when the author was a graduate student at 
the University of Texas at Austin. The author wishes to thank Craig Wheeler for helpful 
comments, John Scalo and Nairn Baliber for encouraging him over the years to write this 
up, and Robert Buchler for giving him the time and financial support necessary to do so. 
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Discussion 
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